Renormalization group dependence of the QCD coupling 
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The general relation between the standard expansion coefficients and the beta function for the 
QCD coupling is exactly derived in a mathematically strict way. It is accordingly found that an 
infinite number of logarithmic terms are lost in the standard expansion with a finite order, and 
these lost terms can be given in a closed form. Numerical calculations, by a new matching-invariant 
coupling with the corresponding beta function to four-loop level, show that the new expansion 
converges much faster. 
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It is of crucial importance to consider the renormaliza- 
tion group (RG) scale dependence of the strong coupling, 
in order to have full consistency in QCD and its applica- 
tions PIQ- As is well known, the QCD running coupling 
a = otsjix = g 2 /(4/7T 2 ) satisfies the RG equation 



da 
du 



(1) 



where the (3 function was calculated to one-loop level 
in QCD more than thirty years ago 0, to two-loop in 
Ref. 13], to three- loop in Ref. p|, and to four- loop in 
Ref. |6J , in the minimal subtraction scheme pj • It can be 
expressed as polynomials of the number of flavors Nf, i.e., 
0i(Nf) = Ylj>o PhjNf. The color factor (3ij is presently 
available to 4- loop level, i.e., 
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with #5,0 = 891C 3 /32 + 149753/768 w 
228.4606573, #5,1 = -1627C 3 /864 - 1078361/20736 « 
-54.26788763, #5,2 = 809C 3 /1296 + 50065/20736 « 
3.164758128, where C, is the Riemann zeta function, 
and C2 = tt 2 /6, Cs ~ 1.202056903, C4 = tt 4 /90, 
£5 « 1.036927755. In this paper, all color factors are 
given for 7V C = 3. Comparing the beta expressions here 
with those in Ref. 6], one would find a difference by a 
factor of 2 2i+1 . 

In practical applications, it is convenient to have an 
explicit expression of the a as a function of the renormal- 
ization point u. The standard approach is to expand it 
to a series of L = 1/ ln(u 2 /A 2 ), where A is the QCD scale 
parameter. However, how the expansion coefficients are 
connected to the beta function is not generally known, 
though one can find the relation to order 3 in Ref. @, 
and to order 4 in Ref. In this letter, the general re- 
lation between the expansion coefficients and the beta 
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function are provided. It is accordingly found that an in- 
finite number of terms like In-' L are lost in the standard 
expansion with a finite order, and these lost terms can be 
given in a closed form. Numerical calculations, by a new 
matching-invariant coupling with the corresponding beta 
function to four-loop level, show that the new expansion, 
with the lost terms included, converges much faster. 
To solve Eq. QJ, let's define 
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With the series expression for (3(a) in Eq. (|TJ, one can 
easily get an explicit expression 
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which indicates that 0(a) is analytic at a — 0, and can 
thus be expanded to a Taylor series as 



(3(a) = J2^ 
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The expansion coefficients $k can be obtained from the 
normal mathematical formula 



^ = h. ^ 



(6) 



Another easy way to obtain these coefficients is to use 
the recursive relation 
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rameter, then Eq. JIJ becomes ,, 



with the obvious initial condition (3q — P2/P1 — (3\/(3q. 
Let L = l/ln(it/A), where A is a dimensional pa- 

-L 2 ^ = (3(a), or 

■^r + ^ + da. Integrating 

this equation gives j — C = + In a+ Jj- fl(x)dx, 
or 
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where C is the constant of integration. 
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Let a = jj-Y(L), then 



— + L* \nY + ^0 k L* k+1 Y k = 1 - CL* - L* In L, (9) 

fc=i 

where L* = (0x/0 o 2 )L, C = {/%/Pi)C - ln/3 , and 
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The graved beta function /3fc can be easily obtained from 
the expression for the acute beta function k in Eq. © 
or (0, and here are the results: 
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where Z = * - 1, flf? = A., flfj = + 

S i ( s = A 8+3 +2, and A = ft/A (i = 0, 1, 2, . . .). 
Now suppose we have a solution of the form 



i=0 .7=0 



(14) 



then 
1 1 
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where the square cup operator, |_|, has been defined in the 
appendix. Substituting these expressions into Eq. J^J, 
we can obtain all fij by comparing the corresponding 
coefficients of u* l \n° u. For — (0,0), (1,0), and 

(1, 1), we have / , = 1, /i,o = C, /i,i = 1. For % > 2 

and j — i, we have X)fc=i( — l) fc Ui o = ®> which gives 
= 1. And for i > 2 and j = i — 1, we get 



(18) 
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whose solution is fi.i-i — iC + X);=i (V — 1) 
For j > 2 and < j < i — 2, we have 
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Please note, there are only terms of fi><i,j'<j on the right 
hand side of this equation. Therefore, it is a recursive 
relation. Here are the solution to order 5: 

h,o = C 2 + C + u (20) 

/ 3 ,o = C 3 + ^ 2 + (3/3i + l)C + /3i+/3 2 , (21) 

fs.i = 3C 2 + 5C + 3/3i + l, (22) 
/ 4 ,o = C 4 + (13/3)C 3 + (6& + 9/2)C 2 

+ (l + 7A 1 +4/3 2 )C + 2/? 1 2 + £ 3 =1 /3 l , (23) 
/ 4 ,i = 4C 3 + 13C 2 + (12/3i+9)C 

+7Ai+4/3 2 + l, (24) 

U,2 = 6C 2 + 13C + 60! + 9/2, (25) 



These correspond to the standard form in Ref. [8j at 
order 3, and agrees to that in Ref. |9j at order 4, i.e., 
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To give a general representation for the expansion co- 
efficients fij, introduce a set of new functions 0i by the 
recursive relation 
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where K = (— l) fc_1 (l + l/A;). From the initial conditions 
(3 = 1 and fli = C, one can easily get all j3i from Eq. 
(127(1 . For C = 0, for example, we have 

Po = 1, h = 0, h = &, h = Pi + 02, (28) 

3 4 



i=l i=l 

h = sLift + y^ + 3/3 2 2 + lW 2 + 5/33, (30) 

Please note, even when one sets all /3i>2 to zero, no Pk>2 
will be zero. This is one of the most important reason 
for us to know the general relation between fcj and the 
beta function. 

It is found that /y satisfies 
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where |Jo ^ are * ne Taylor coefficients of ln fe (l + x), 
i.e., ln fc (l+x) = E^o(US*0 X ' = (E^*^)* = 
(y^Li jx l ) k ■ Accordingly we have <f>o = 0, $;>i = 
l/l, Uo ^ = 4o, Uo $ o = 4,0, Uo ®l = ®h an d 

u*«= n e 
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One can naturally consider to prove Eq. I|31|l by math- 
ematical induction. We have numerically checked it to 
order 100. Now we directly used it to give an expression 
for fij, i.e., 
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where P^j. can be obtained by repeatedly using Eq. (|31|l . 
For fe = i — j and fe = i — j — 1, it is easy to get 
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Generally for k = i — j — I or I = i — j — k, careful 
derivations give 

P}$ = Q(i,j,l)Q(i- j,l,Q) 

l-l ( ( r l-l-ESN.N 

+E m E [e(u,-x) 

r=l I p= i a p =l / 

x Q (i - j - SgZ^Sg, s r , Z - SJ =1 s 9 ) 
xe(i-i-E^i«„I-E^i« g ,0) 

x II ^ r ~ ~ s ?=i s 9' vk) p (35) 



where Q is a function of three non-negative integers and 
defined by 
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matrix of order i + 1. The specially simple elements 
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In the traditional minimum subtraction scheme (MS), 
the strong coupling a(u) as a function of the re- 
normalization point u is not continuous at the quark 
masses. Let's derive a matching-invariant coupling by 
absorbing loop effects into the MS definition and give 
the corresponding beta function to four-loop level. 

Suppose the new coupling a' is connected to the orig- 
inal coupling a by 



i=0 



Then, using the matching condition a = Ej=o w j 
with the matching coefficients 0, 
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where an overhead check means decreasing Nf by one 
flavor to the corresponding (N{ — l)-flavor effective the- 
ory. Accordingly, comparing the coefficients of a in the 
equality a' — a' yields 
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Assume a t = Y?j=o a i,j^{ \ tncn °>k = Z)j=o a k,j{N{-iy 
Substitution into Eq. 142|) then gives 
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whose solution is 

a Q = &0,0) a l = a l : 0, 0-2 = ^2,0 

Y7037 82043 



^a , N u (44) 



(':>, 



^3,0 " 

2633 



V 1536 27648 



11 



62208 



f > 



N{ 
(45) 



To definitely fix the new coupling, one needs to choose 
a, o- The simplest choice would be a i — 5 i . With this 
convention, one has 
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where a 3 ,i = 7037/1536 + 82043C 3 /27648 « 8.148377983. 
Then the new matching-invariant coupling is 
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The renormalization equation for a' is 
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The primed beta function ft' can be obtained as such. 
Operating with u-^ on both sides of Eq. (|38|l . applying 
Eqs. and (Q), and then comparing coefficients will 

give XTfc=o P'k Uo +2 - (fc + l)afeft_fc = 0, namely, 
ft are given by the recursive relation 
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On application of Eqs. 1)46(1 and 1)47(1 . one immediately 
has the following explicit expressions for the new beta 
function: 
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180 
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with ft 

ftf 5 = 2118091/82944 + (2603291/55296)C 3 ~ 

82.12826576, (3f ] = -4027/1296 - (194353/82944)C 3 « 
-5.923892810. 

It should be mentioned that a different expression for 
ft was previously given in Ref. 11]. The difference is 



caused by the fact that a wrong value for C 2 was quoted 
there Gj. 

As an application of the general relation between fij 
and the beta function, one can develop another expansion 
which converges much faster. For this one can observe, 
more carefully, the standard expansion 
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Representing the terms in this expansion with the corre- 
sponding coefficients fij, all the terms can be arranged 
in a matrix as 
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The standard expansion corresponds to summing the 
terms row by row. When one takes the expansion to a 
finite order, i.e., replacing the oo in Eq. ((551) with a pos- 
itive integer N, as has been done in the usual way, then 
all the terms like fjj (N < j < oo) on the diagonal and 
fj+ij on the next to diagonal are missed, although these 
terms are all known and have nothing to do with beta 
functions. Generally, the terms fj+k,j for < j < oo on 
the kth next to diagonal involves only ft<z<fc_ 2 . But all 
the terms fj+kj with j > N are lost, though no such 
terms are zero even when one sets all ft>2 to zero. 

To include the contribution from the terms just men- 
tioned, we can consider to sum over diagonals, which can 
be achieved by taking i = j + k in Eq. ((55(1 . i.e., 
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where the expressions for fj+k,j can be obtained from 
Eq. J23J): 
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From these expressions, we can give compact form to Ik- 
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where x = 1 + (fli/Pl) lnln(u/A)/ln(u/A), and X = 
L*/{0xx) = l/tflg ln(u/A) + ft lnln(u/A)]. 




o.o r .. .. i . — ■ — . — 

1 10 100 

u (GeV) 



FIG. 1: The QCD coupling as functions of the t'Hooft unit 
of mass. The dot, dash, solid, and dot-dash lines are, respec- 
tively, for the order TV from 1 to 4 with the scale A indicated in 
Tab. 1. The full dots are the experimental data of, in increas- 
ing order of A, r width, T decay, deep inelastic scattering, 
e + e~ event shapes at 22 GeV from the JADE data, shapes 
at TRISTAN at 58 GeV, Z width, and e + e" event shapes at 
135 and 189 GeV. 
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FIG. 2: The relative difference between the new expansion 
in Eq. I|63fl and the conventional expansion in Eq. I|55fl . 

InEq. there are an infinite number of logarithmic 
terms like hr' x which are included in Ik , even when one 
takes the expansion to a finite order, say a = X^fcLo^- 
The first several IkS have been worked out to a closed 
form in Eqs. (ftT0")l - (frT2*)> . The procedure is, to some ex- 
tent, similar to that in Ref. 13] for the matching func- 
tion. Here we can, in fact, give the closed form for Ik to 
arbitrary order: 

oo k 

a = J2 P^i Xk+1 E^ 1 )'^ ln ' x - ( 63 ) 
fc=0 z=o 



In Eqs. I|63(l and (|55|l . there are two arbitrary con- 
stants: A and C. Because the renormalization group 
equation 1)49(1 or Q is of the first order, only one of 
them is independent. So we can arbitrarily take one of 
them, while the other is determined by giving an ini- 
tial condition. It nearly becomes standard, nowadays, 
to take C = [Tij . which makes expressions somewhat 
simpler, and a(mz) = 0.1187/7T, where Mz = 91.1876 
GeV is the mass of Z bosons. Setting C = requires 
distinct A for different effective flavor regimes, and we 
use A 6 , A 5 , A 4 , and A 3 for u > m t , mi, < u < m f , 
m c < u < rrifc, and m s < u < m c , respectively, where 
the relevant quark masses are taken, in the present cal- 
culations, to be ra t — 175 GeV, m& = 4.2 GeV, m c = 1.2 
GeV, and m s = 100 MeV. 

In Fig. ^ the coupling is shown as a function of the 
renormalization point, calculated from Eq. (|63|l with the 
infinity replaced by N, and the order N from 1 to 4. The 
same calculation has also been performed from the con- 
ventional expansion in Eq. J5SJ. The relative difference 
between the results from Eq. (|(j3l) and Eq. (|55|l is shown 
in Fig. |21 It can be seen that, with decreasing u, the 
difference becomes more and more significant. 

TABLE I: QCD renormalization group scale parameter A for 
the order from 1 to 4. For each Ai, the left column is for 
the new expansion in Eq. (1631 and the right column is for the 
conventional expansion in Eq. 15511 . 
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To compare the convergence speed of Eqs. (I63f) and 
(|55|l . all the A^ (i — 3-6) are listed in Tab. D There 
are two columns corresponding to each A i; the left col- 
umn is for Eq. (|63|l while the right column is for Eq. 

. It is obvious that the new expansion ijfj3l) converges 
much faster than the original expansion (|55|l . Even at 
the leading-order (N — 1), the corresponding hi for Eq. 
(iBl^l has nearly approached to its value at order 4. So 
in practical applications, it should be very accurate to 
calculate the coupling simply by 



^ln(u/A) + /3ilnln(u/A)' v ' 

In summary, the general relation between the standard 
expansion coefficients and the beta function is carefully 
derived. A matching-invariant coupling is given with the 
corresponding beta function to four-loop level. A new 
expansion for the coupling is then deduced, which is in 
principle more accurate than the conventional expansion 
due to the inclusion of an infinite number of logarithmic 
terms in a closed form. 
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APPENDIX A: SQUARE CUP OPERATOR 



The square cup operator, |_|^, is denned so that it 
meets 



E°« 



X |J"< (A1) 
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Obviously one has |J m a { = 5 ifi , \J m a l<m = 
0, \_^ m fti> m = aj. And for k > 2, it is also not difficult 
to give a general explicit expression 
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where 
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Here are special examples: 



where 



^jpift * ->1 . (A3) 
otherwise 
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The meaning of <r£ is similar to this. Here are several 
special simple cases: 



[Jao = £fc,o, LJa*<fcm=0, |Ja fem = a^, (A4) 



The two-dimensional extension of the square cup op- 
erator, |_| m n , is defined by 



oo i \ oo i / k \ 

EE EE U/<. ( A5 ) 

,i=mj=n J !=0 3=0 \m,n / 
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